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Abstract: We link the FDR treatment of ultraviolet (UV) divergences to dimensional 
regularization up to two loops in QCD. This allows us to derive the one-loop and two-loop 
coupling constant and quark mass shifts necessary to translate infrared hnite quantities 
computed in FDR to the MS renormalization scheme. As a by-product of our analysis, 
we solve a problem analogous to the breakdown of unitarity in the Four Dimensional 
Helicity (FDH) method beyond one loop. A fix to FDH is then presented that preserves 
the renormalizability properties of QCD without introducing evanescent quantities. 
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1. Introduction 

When using customary approaches [1-5] to handle the ultraviolet (UV) problem in quan¬ 
tum field theory (QFT), intermediate steps are necessary to extract physical answers from 
loop calculations. In particular, the Lagrangian C of the theory is modified by adding UV 
counterterms (CTs). They absorb the divergences generated by the high-frequency part of 
the loop integrals and as one moves up the orders in the loop expansion one must include 
lower loop CT calculations in order to be consistent. When carrying out this renormaliza¬ 
tion program UV infinities are usually regulated via dimensional regularization [4] (DReg) 
and renormalized quantities are defined by specifying what is subtracted from the bare 
ones. A particularly convenient subtraction scheme is MS, in which only UV poles and 
universal constants are dropped. 
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The FDR^ approach of [ 6 ] deals with UV infinities in a different way. A new kind of 
loop integration^ is introduced that coincides with Riemann integration in UV finite cases, 
but produces a finite and regulator free answer also when acting on divergent integrands.^ 
In this way no CTs need to be incorporated into T: they are traded for a change in 
the definition of the loop integration. Moreover, FDR directly generates renormalized 
amplitudes since it is independent of any UV cutoff. 

The main aim of this paper is to construct the one- and two-loop transition rules 
between FDR and MS in the framework of QCD. As a perturbative treatment of a renor- 
malizable QFT is unique up to a renormalization scheme dependence, we have achieved 
this task by studying off-shell^ Uloop QCD correlation functions and by searching for 
a DReg renormalization scheme which reproduces the FDR answer when i = 1 and i = 2, 
i.e. 


^renormalized — ^'Lre + (Aloop-CTs) (1-loop-CTs), (1.1a) 

^renormalized ^FDR‘ 

We dub such a scheme Djfeg^ and its renormalization constants These can then 

be used to extract the coupling constant and quark mass shifts that relate FDR to MS. 
We emphasize that in a typical FDR calculation there are no since one directly 

computes eq. (1.1b). We introduce them in the context of this paper only because we want 
to work out the correspondence of FDR with a canonical renormalization approach based 
on counterterms. 

The second insight of our work concerns the nature of self contractions of Lorentz 
indices and 7 -matrices explicitly appearing in denote by Ug. When studying 

eq. (1.1) one needs to reconcile the value Ug = n dictated by DReg in (1.1a) with Ug = 4 used 
in (1.1b). This has to be done without spoiling the renormalizability of QCD. Following 
FDR as a guidance, we have been able to establish the correct —)• 4 limit of DReg, 

compatible with a local subtraction parametrized in terms of the A mismatch 

between the structure of G^are the CTs has been recently recognized [7-9] to cause the 
failure of the renormalization program beyond one-loop in the original formulation of the 
Four Dimensional Helicity scheme (FDH) [10,11]. Since we observe that terms restoring 
the right cancellations are produced by FDR, we argue that the FDR approach provides a 
natural solution to this problem. In contrast with methods that require the introduction 
of e-scalars [ 12 ], in our case no new fields nor evanescent couplings need to be added to the 
QCD Lagrangian. Moreover, a shift linking the MS quark mass to a fixed version of FDH 
can be worked out. So far we have tested this approach for off-shell correlation functions. 
What happens in the presence of IR divergences needs to be further investigated. 

^Acronym of Four Dimensional Regularization/Renormalization. 

^Called FDR integration. 

^Abusing a bit the language, we dub divergent (convergent) integrands those which would generate UV 
divergent(convergent) integrals upon normal, four-dimensional integration. 

"^Working off-shell allows us to deal with IR convergent integrals. A detailed two-loop study of FDR in 
the presence of IR divergences is equally important, but outside our focus. 
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The outline of the paper is as follows. In the next section we review FDR. Sections 3 
and 4 describe our computational strategy. Our results are collected in section 5 and 
appendix A. Furthermore, in section 6 we discuss our fix to FDH. Finally, appendices B-D 
contain explicit examples of our algorithm. 


2. FDR 


FDR integration was first introduced in [6] and several examples of one- and two-loop 
computations can be found in [13-17]. Here we briefly review its definition, which we need 
to make contact with our calculation. 

Consider an £-loop integrand J{qi,... ,q£) depending on £ loop momenta qi,... ,q(. 
The multi-loop FDR integration over J is defined as^ 


/ 


[d'^qi] ■ ■ ■ [d'^qe] J{qi,...,qe, = lim / d'^qi 


' d'^dedpigi, ■ ■ ■ ,qe, ( 2 . 1 ) 


where ..., qe, is the UV finite part of J{qi,... ,q£, fj?), specified below, and fi a 

vanishing mass required to extract Jp from J. J(gi,..., ®, and Jf(( 7 i, ..., q^, are 
read off from J{qi,... ,qi) according to the following rules: 


i) Squares of integration momenta appearing both in the denominators of J(gi,... ,qi) 
and in contractions generated by Feynman rules are shifted by 

q^ ^ qf - ^ ql ( 2 . 2 ) 

This replacement is called global prescription] 

ii) A splitting 


J{qi,...,qe,H^) = [JwF{qi,...,qe,n^)] + JF{qi,...,qe,ia^) (2.3) 


is performed in such a way that UV divergences are entirely parametrized in terms 
of divergent integrands® contained in [Jinf], that solely depend on p?; 

hi) The global prescription of eq. (2.2) should be made compatible with a key property 
of multi-loop calculus:^ 

in an i-loop diagram, one should be able to calculate a sub-diagram, 
insert the integrated form into the full diagram and get the same answer. 

We dub this sub-integration consistency. 

®FDR integration is denoted by the symbol [d^Qi]. 

®By convention, we write divergent integrands between square brackets and call them FDR vacua, or 
simply vacua. Examples of the extraction of FDR vacua from loop integrands are given in appendix B. 

’^This last requirement turns out to be the mechanism that enforces the renormalizability of the Daeg^ 
scheme. It is discussed at length in section 4.3. 
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Finally, after lim^_^o is taken, In^ —)■ In^^ is understood in the r.h.s. of eq. (2.1), where 
//r is an arbitrary renormalization scale. 

This definition preserves shift invariance 

J [d^qi] ... [d'^qe] J{qi,.. .,qe,d^‘^) = J [d^Qi] ■ ■ ■ [d^Qe] J{qi + Pi, ■ ■ ■,qe+Pe,P^), (2-5) 

and the possibility of canceling numerators and denominators 

fid*,!,] ■ ■ ■ id%] = Jid^d,] ■ ■ ■ 

which are properties needed to retain the symmetries of C. From eqs. (2.5) and (2.6) it 
follows that algebraic manipulations in FDR integrands are allowed as if they where conver¬ 
gent ones. This authorizes one to reduce complicated multi-loop integrals to a limited set of 
Master Integrals (MI) by using four-dimensional tensor decomposition [14] or integration- 
by-parts identities [17]. In other words, the definition in eq. (2.1) can be applied just at 
the end of the calculation, when the actual value of the Mis is needed. 

An important subtlety implied by eq. (2.6) is that the needed cancellation works only 
if integrands involving explicit powers of in the numerator are also split via eq. (2.3) 
as if = q^, where q^ is the momentum squared which generates jjf. As a consequence, 
although only one kind of jj? exists, one has to keep track of its origin when it appears in 
the numerator of J{qi,... ,q(, iJp). For this we use the notation which understands 
the same splitting required for qf. FDR integrals with powers of in the numerator are 
called “extra integrals”. Their computation is elementary. One- and two-loop examples 
can be found in [6,14] and appendix D. 


3. FDR and DReg on the same footing 


Eq. (1.1) defines the Djfeg^ scheme we are looking for, i.e. the DReg scheme reproducing 
the FDR correlation functions. As the r.h.s. of (1.1a) is computed in DReg while the r.h.s. 
of (1.1b) in FDR, we need a common framework which accommodates both approaches. 
Here we illustrate this framework at two loops, but the same considerations apply at any 
loop order. 

Our starting point is eq. (2.3) when J{qi, ( 72 ) is the sum of all integrands contributing 
to Gbare- V J{qi,Q 2 ) IS free of IR infinities we can rewrite the DReg integration over it as® 


/ 


d'^qid'^q2J{qi,q2)= / d'^qid'^q2 lim J(gi, ^ 2 , M ) = 1™ / d'^qid'^q2 J{qi,q2, P ) 

J J 


= liin / dr'qid^q 2 [JiNF{qi,q 2 , + hni / d‘^qid‘^q 2 JFiqi,q 2 , 


fjj —^0 


11 —^0 ^ 


(3.1) 


where the absence of IR divergences authorizes us to extract lim^_^o out of the integral. 
The integral over J-p{qi,q 2 , is and through eq. (3.1) we have isolated it within 

®n is the number of dimensions defined as n = 4 — 2e. 
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^bare’ Scheme can be determined by solely looking at the pieces which 

differ between eq. (1.1a) and (1.1b), namely 

[«^inf(<?i, <? 2 ,1^^)] , (1-loop-CTs) and (2-loop-CTs). (3.2) 

This defines our strategy. For instance the (1-loop-CTs) are known at two loops, so that 
Djfeg^ is defined by choosing 

(2-loop-CTs) =-(1-loop-CTs) - lim /[</iNF(gi, 92 ,Z^^)] , (3.3) 

which sets the finite part of G^^re FDR value. 

An important condition should be fulfilled by eq. (3.3). It should be local and propor¬ 
tional to the the Born correlation function G®: 

(2-loop-CTs) = Const x G°, (3.4) 

which is the correct form to be re-absorbed into We dub eq. (3.4) our renor- 

malizability condition. To understand its consequences it is convenient to split [Jinf] as 
follows; 

[Jinf{qi,Q2, C^)] = [<^ 11 )^( 91 ) 92 ,/^^)] + [<^ 1 ^( 91 , 92 ,/ i ^)] • ( 3 . 5 ) 

[J^j^f] is a factorizable contribution in which only one sub-integration is UV divergent. 
It is called the sub-vacuum and depends on the kinematic scales entering the UV finite 
sub-integration. The second piece is independent of kinematics and it is called the global 
vacuum.^ Eqs. (3.3) and (3.4) require that, upon integration, kinematic dependent terms 
in [e/f^p] should cancel one-loop counterterms. In section 4 we will describe the subtleties 
of this cancellation. 


4. The calculation 


As outlined in the previous section the core of our strategy is calculating FDR vacua in 
DReg instead of throwing them away. As a first step we write down^® the twelve integrands 
J'-'i (i = 1,..., 6 and / = 1,2) corresponding to the one-particle irreducible one- and two- 
loop QCD correlation functions drawn in figure 1, excluding CT diagrams. Then, we split 
them as in eq. (2.3) and compute the integrals over their vacuum parts. 


// = lim 


d^q 


rG 


''^inf(9, 


and 


If = lim 


dl'qi dP'q 2 


rG\ 


''^inf(9i, 92 , h- ) , (4.1) 


in the symmetric off-shell kinematic point 


kl = kl = kl = M^. 


(4.2) 


With this choice all integrals are made free of IR divergences, so eq. (3.1) applies. Further¬ 
more, the complexity of the calculation is reduced to a one-scale problem. An important 

^Examples of global vacua and sub-vacua are given in the first line of eq. (B.15) and in eq. (B.18), 
respectively. 

^°We work in the Feynman-’t Hooft gauge using QGRAF [18] for generating the diagrams and FORM [19] 
for extracting their vacuum part. 
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Figure 1: One-particle irreducible one-loop {£ = 1) and two-loop {i = 2) QCD Green’s functions 
used in our calculation. The gray blobs denote the sum of all possible Feynman diagrams computed 
in the Feynman-’t Hooft gauge. Diagrams with counterterms are not included when the are 
used to compute the r.h.s. of eq. (1.1b). 


point concerns explicit contractions appearing in the numerator of 



7“7a = = ns- 


(4.3) 


As 1} and 1“^ are regulated in DReg, rig = n should be used in eq. (4.1). 

The knowledge of the allows us to parametrize the FDR subtraction in terms of 
renormalization constants However, not all of them are independent, as they are 

related by QCD Slavnov-Taylor (S-T) identities.In the following we consider, in turn, 
the 1 = 1 and 1 = 2 contributions to the 


4.1 The one-loop case 

At one loop the calculation of is simple. Schematically: 

• The global prescription of eq. (2.2) ensures that reducible g^s in the numerator sim¬ 
plify with denominators; 

• The integrands are then split as in eq. (2.3) by using the propagator identity 

(4.4) 

(4.5) 

^^See section 5. 


1 


1 


2(g ■ k) + 


+ 2{q ■ k) + k"^ q^ q^{q^ + 2{q ■ k) + k"^) 
which allows one to express in terms of tensors of the kind 


d^q q<^2r 

^ ’ 


(4-b2r-2s > 0); 
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Finally, by virtue of lim^_>.o, polynomially divergent integrals vanish, so any tensor 
structure is reduced to a fundamental scalar 

f d^q q°‘^ ■ ■ ■ 2 


with 




^-2^ (g 2 - ^ 2 )r +2 (2r + 2)!! 

d'^q 




1 






1 


A =- 7 £; - In TT, 

e 


(4.6) 


(4.7) 

(4.8) 


and where ^"i'"" 2 r jg completely symmetric and only made of products of metric 
tensors. 


Eventually, the logarithm in eq. (4.7) has to be combined with the one-loop analogue of 
the Jf term in eq. (3.1) to compensate the ^ dependence of the finite part.^2 As a result, 
the UV divergent part of Jl, Inf(/j^), is fully proportional to^^ 

Vo=Viiin)l^o = i7r^A. (4.9) 


Now we study Fin(/1), namely the UV convergent contribution to //. It determines the 
finite part of the one-loop subtraction term in eq. (1.1a) and appears due to the explicit 


presence of Ug in J^p{q,IJ- 


. In fact, rewriting 
= 4 — 2Ae, 


(4.10) 


we see that when Ae multiplies UV single poles, constants are generated that are fully 
subtracted in Difeg^ but contribute to the finite part of Gj in MS. Thus, the requirement 
of eq. (1.1b) causes a deviation from the minimality of the zj' proportional to (1 — A), 
such that A = 0 (1) in Djfeg^ (MS). As one finds that both Inf(/1) and Fin(/1) factorize 
the Born, one computes 


^1,FDR^0 ^ _jl 


(4.11) 


In summary, at one loop it is possible to perform a calculation in DReg and consis¬ 
tently renormalize the result to reproduce the FDR answer. In other words, FDR can be 
interpreted as a particular renormalization scheme of DReg, i.e. the Djfeg^ scheme we are 
looking for. 


4.2 The two-loop case 

Here we illustrate the calculation for massless QCD. When dealing with a non vanishing 
quark mass niq the formulae complicate a bit, but the reasoning remains unchanged. For 
simplicity, we start from cases^^ in which we see that the sub-integration consistency (2.4) 
does not play any role, and postpone to section 4.3 the study of more complicated situations. 

After scalarization by means of tensor reduction and integration-by-parts, one finds 
two types of contributions to if: 

^^This is why FDR integration is defined with the replacement In/r —> In/iR. 

^^The notation means neglecting terms of 0(e). 

^"^Specified later on. 
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An integral over the global vacuum of the kind 


GVi(fi) = lim 

fi —^0 ^ 




Fii{n,ns 


1 


-4 -4 


+ Fi2{n,ns 


1 


v<ii<ihi2 


(4.12) 


where Fij{n,ns) are rational functions and qu = Qi + Q 2 ] 
An integral over the factorizable sub-vacuum^^ of the form 


SVii^i) = lim / 
/i-j-O J 

in which Ji{qi,fj,) is UV convergent. 
By judiciously using the identity 


d'-^qi <Fq2 

.-2e ,,-2e 


/^R /^R 1^2 




(4.14) 


1 


qf + 2{qi ■ k) + k‘^ 


1 

¥ 


qf + 2{qi ■ k) 


k?{ql + 2{qi ■ k) + k"^) 


(4.15) 


Ji{qi, fj,) can be split into a piece which develops a In^u^ upon integration and a term where 
H can be set to zero at the integrand level 


¥ = Jaiiqi, ¥ + Jbiiqi)- 

With our special kinematics we find 

Ai{n,ns) 


(4.16) 


Jai (^ 1 1 ¥ — 


Jbi{qi) — 


qf 


Ai{n,ns) 

qKqi + kiY 


+ M" 


Bi{n,ns) 


qiiqi - kiyiqi + k2y' 


(4.17) 


where Ai{n,ns), A'-{n,ns) and Bi{n,ns) are rational functions. Notice that the UV finite¬ 
ness of Ji{qi,fi) ensures that the the pole parts of the integrals over Jai{qi,¥ and Jbiiqi) 
cancel each other. Furthermore, Bi{n, Ug) = 0 in the case of two-point correlation functions. 
An explicit example of the procedure yielding eqs. (4.12)-(4.17) is reported in appendix B. 
The additional complexity at two loops is that, unlike 


the integral 


SVaiifi) = lim 

l-L —^0 ^ 


SVbi = lini 

/i—^0 ^ 


d^qi dAq2 1 

[qt 


dAqi dFq2 

- 2 t ,- 2 e 


Mr L^2J 


Jaiiqij l-d), 


Jbiiqi') 


(4.18) 


(4.19) 


®This sub-vacuum is generated by the application of the identities (with i A j) 


1 


(lij - 2(?j • 90) 


1 


1-t 


-I- 2{qi ■ qj) 


1 


^ _ 9j -t 2(gi ■ qj) 


(4.13) 


needed to disentangle the sub-divergences of if from the finite part. 



depends on the kinematics. Thus, it generates logarithms of physical scales that cannot 
be absorbed into the Therefore eq. (3.3) tells us that, in order to establish the 

connection between FDR and the standard renormalization approach, we must demonstrate 
that such non-local terms are compensated by the sum of all diagrams containing insertions 
of one-loop counterterms.This contribution is dubbed CTi in the following. To achieve 
this cancellation, we have to recast eq. (4.19) into a form suitable to be combined with the 
CTi. Observing that^^ 


/ rl^n 

^J,((?i,/i) + 0(e) 

/^R 


leads us to consider 


(4.20) 


SVM = Volhn f ^J,,(gi,^), 

/ d^ a 

(4.21) 

/^R 

instead of eqs. (4.18) and (4.19). Now has the same structure of a counterterm 
diagram, as it factorizes Vq. An explicit calculation of the CTi gives 


CTi + SVli = 0 for i = l, 2,4,5, 


(4.22) 


and we obtain, for correlation functions without external fermions, the following two-loop 
renormalization constants 


,72,FDR 






(4.23) 


in accordance with the renormalizability condition of eq. (3.4). As in eq. (4.9), the choice 
of the point fi = has the effect of removing the dependence on ij, from the last integral 
in eq. (3.1). 

As in the one-loop case, Ug in eqs. (4.12) and (4.17) generates hnite terms proportional 
to (1 —A) in the In addition, finite contributions are created by two-loop integrands 

that do not multiply powers of Ae, e.g. 


fjn) 


(4.24) 


To denote their origin we multiply by a parameter 5 any finite combination 



fin) 

-4 -2 -2 

— Pole Part < 


fin) 

-4 -2 -2 

/ 

191^2^12] 

1 

J 

lqtQm2\ 


(4.25) 


which is fully subtracted in but contributes to the finite part of G? in MS, so that 

5 = 1 (0) in (MS). 


16 1 FDR 

For consistency with the form of the Z^ ’ the one-loop insertions have to be computed with A = 0. 
is defined in eq. (4.9). 
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4.3 “Extra”-extra integrals and snb-prescription 

In the case of QCD Green’s functions with external fermions eq. (4.22) does not hold true^® 
and leads to results incompatible with eq. (3.4). In this section, we show that the FDR sub¬ 
integration consistency (2.4) requires the introduction in the finite part of the correlation 
functions of “extra”-extra integrals {EEIs) with the same structure of SV/, namely 

EEIi[^) = EEIai{n) + EEIhi. (4.26) 

They are just what is needed to restore the renormalizability condition, i.e.^® 

CT, + - EEIaiifi) - EEhi = -EEIM (4.27) 

and 

^2,FDR^0 ^ + EEIaM {i = 3, 6), (4.28) 

with EEIai independent of physical scales. 

The FDR origin of the EEIi is the need of introducing a sub-prescription to cure the 
mismatch between the global prescription of eq. (2.2) and the consistency condition (2.4). 
In fact, although essential to preserve gauge invariance at two loops, the shift 

Qi ^4- Ai (4-29) 

may clash with the analogous replacement required to ensure (2.4) at the level of divergent 
one-loop sub-diagrams. This is better explained with an example. Consider the two-loop 
diagram^*^ of figure 2. Before applying eq. (4.29) its numerator reads 



Figure 2: The two-loop diagram contributing the Nf corrections to G§. N{qi,q 2 ) is given in the 
text dropping irrelevant constants. The replacement qf —> q^ is performed only in the denominators. 

Niqi,q2) = 47a(^i -b h)lp • ^ 12 ) + 9?292 + (I 12 Q 2 ) ■ (4-30) 

The sub-prescription is defined as the effect of eq. (4.29) on N{qi,q 2 ) from the point of 
view of the divergent sub-diagram pulled out from the rest of the diagram, as in hgure 3. 

From the perspective of the q 2 sub-integration one has to distinguish internal and 
external parts. We denote the separation of the parts external to the sub-diagram by 
placing hats on the Lorentz indices. Algebraically the hats do not make any difference 
other than to denote the fact that they are of an origin which is external to the sub¬ 
diagram, so all standard identities apply, for example 

_ = 7a7" = 4, 7«92“ = h- (4-31) 

^^Interestingly, it applies when computing if and 7| with Us = 4, which is the value prescribed by FDH, 
as discusses in section 6. Nevertheless, we emphasize that we must set Us to n because our strategy is to 
reproduce the FDR result by regularizing FDR vacua in DReg, which dictates ris = n. 

'^^When calculating the vacuum EE Is contribute with a minus sign and should be computed with Us = n. 
■^^This is the only possible contribution to G3 proportional to Nf. 
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We then study how the part^^ of N{qi, ^ 2 ) divergent in q 2 

+ h)Q2 + + h)h, (4.32) 

transforms under eq. (4.29). From the perspective of the sub-prescription the second term 
is inert, since it has an origin external to the sub-diagram. Thus 

N^^Hqi, q2) ^ N^^\q^,q2) - 8(^1 + (4.33) 

On the other hand, the two-loop global prescription dictates that also the second term 
must transform as^^ 


+ h)i2 H2{ii + h)h + 8(^1 -h h)^J‘^\2, 


leading to no change in 


N^‘^\qi,q2). 


(4.34) 

(4.35) 


Therefore, to remove double counting, the EEI to be added to the diagram in figure 2 is 
defined as the result of the sub-prescription minus the outcome of the global prescription, 
i.e. the difference^^ between the r.h.s. of eqs. (4.33) and (4.35) 

+ h) 


EEI = -I 


! J [d^gi][d^g2]- 


(4.36) 


IiQhi 2 {Qi +kf + 2{qi ■ ki))' 

Here, we mark the fP \2 with a hat because it is acting only on the q 2 sub-integral. One 
finds^"^ 


EEI{^jO = -iTT^} 




1 


= Ho 


4 — n, 


q^{q^ + k‘l + 2{q ■ ki)) 
(Eq 1 


— EE If) EEIa{iJ,ji) 


— lim 


iaj^‘^^q‘^{q‘^ + kl + 2{q-ki)) J q 


f (Eq 

1 

1 

/ -2e 




9 

II 


(4.37) 


where the last representation is suitable to be used to prove eq. (4.27). A further example 
of computation of an EEI via sub-prescription is given in appendix C. 

Finally, it is important to realize that 
EE Is also arise in the intermediate steps 
of the calculation of the QCD correlation 
functions without external fermions dis¬ 
cussed in section 4.2. From eq. (4.22) and 
the universality^® of the coupling constant 
extracted from all of the three-point ver¬ 
tices in figure 1 we infer that one must find 
zero when summing up the EE Is from all 
contributing diagrams. Indeed, we have 
explicitly checked that this happens in the 
case of the Nf corrections to the gluon and ghost propagators. 

^^Here this is quadratic because less-than-quadratic terms do not generate extra integrals. 

^^This is obtained by anticommuting until the q 2 meet. 

doing so, we reintroduce the correct one-loop behaviour leading to (2.4). 

^“^See appendix D. 

^®See section 5. 



Figure 3: The same diagram of figure 2 from the 
point of view of the (72 sub-integration. Lorentz 
indices external to the the sub-diagram are given 
a hat. 
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5. Results 


The main result of our calculation is the list of QCD renormalization constants in ap¬ 
pendix A. When setting (A, 5) = (1,0) the known MS formulae [20-22] are recovered. 
The choice (A, 6) = (0,1) corresponds to Djfeg^, i.e. the renormalization scheme which 
reproduces the FDR QCD Green’s functions up to two loops. 

As required by the QCD S-T identity 


ZgGG Zgcc /— — j— — 

- = \/ ZasV ZgG, 

^GG ^cc 
we have verified that the three ratios 

1 f Zggg ^^ 1 f Zgcc '\^ _J_ Zg^\^ 
ZgG V Zgg J ’ Zgg \ Zee J ’ Zgg V Zq,qf ) 


(5.1) 


(5.2) 


produce the same coupling constant renormalization, namely Zog(A,(5) in eq. (A.9). In the 
case of Zag{0, 1) this provides a stringent test on the universality of FDR. 

Using niq 7 ^ 0 in Gg gives the quark mass renormalization constant in eq. (A. 10). 
Again the correct MS result [ 21 ] is reproduced with A = 1 and 5 = 0. As an extra check, 
we have verified that any niq dependence drops in Zq,q,, as should be. 

Eq. (A.9) can be used to infer the coupling constant shift between FDR and MS 


Z„g(l, 0 ) a 


FDR 


ZasiOA) 


= = 1 + 


a 


MS 


a 


MS 


dvr 


Nc 

— + 
3 


a 


MS 


dvr 


-N^ + 8N^f 


+Nf 


Nc -^Cp-f 


2 d 

-^Nc + -Cp 


(5.3) 


In an analogous way eq. (A. 10) produces the quark mass shift 

2 


= 1-Cf( ^ ] + C'f ( ^ ] \ —Nc - -Cf + f(9Nc + —Cf 
MS I dvr / M dvr M 2 d 8 V 3 


rrii: 




(5,4) 


Eqs. (5.3) and (5.d) provide the transition rules from IR finite QCD quantities computed 
in FDR and their analogue in MS. 

Finally, it is well known that one can use the coupling constant shift between two 
schemes at a given order to relate the two beta functions at one order higher. In the 
following we sketch out the two-loop proof of this [ 11 ] and use it to derive the three-loop 
beta function in FDR. Let us consider two different renormalization schemes defined in 
terms of two coupling constants a a and ap related by a shift 


otA 

dvr 


OiB 

dvr 


1 -I- Cl 




(5.5) 
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The two beta functions 


d aA,B {OiA,BY uA,B , i^A,B fO:A,B\ , ^A,B (OtA^B^ , ^ z' 3 ^ 

= "5;;^ = (u^rJ "» +"■ (u^rJ +0 (“-’.«) 


dfi 47r 


(5.6) 


are linked by the chain rule 


d a A ( d a b\ da a 

^ ^djj, Att ^ \d^ Att ) daB 


daA 

PBj —• 

daB 


(5.7) 


Writing both sides in terms of the same a through eq. (5.5) we find the standard result of 
scheme independence up to two loops, i.e. and bf = hf as well as an expression 

for the three-loop beta function in scheme B, using only two-loop information from scheme 

A 

b2 =b2 + cibf + {cl - C 2 ) 60 . (5.8) 

To calculate the three-loop beta function in FDR we use the three-loop MS beta function 
in [23] together with the values from from eq. (5.3). This gives 




= Nil- 


3610 176 


+Nf 


27 

1331 


+ -^N, + f(~CF- 


9 


27 


9 


:Nf 


n 292 o 


7140^0 88 

—N^ + — 




(5.9) 


6. Fixing the UV behaviour of FDH 

Recently it has been observed that the original formulation of the FDH method breaks 
unitarity [7]. This was seen as the standard renormalization program fails to remove all 
UV poles at high enough perturbative orders. During our work, it proved of interest to 
try and investigate this problem from our renormalization scheme perspective. Using the 
technology developed in the previous sections we can easily perform an analysis of the 
validity of the interpretation of FDH as a perturbative description of QFT. 

We begin by setting up the problem as in eqs. (1.1). We first recall eq. (1.1a) for an 
arbitrarily renormalized, dimensionally regulated correlation function. Next, we equate it 
to correlation functions calculated in the FDH scheme up to two loops, i.e. 

T^fenormalized = + (Uloop-CTs) (1-loop-CTs), (6.1a) 

^^renormalized {7rbare(^.s ^)} ’ ( 6 . 1 b) 

^^renormalized = MS {GLe(^. = 4) + (l-loop-CTs)U.= 4 } . ( 6 . 1 c) 

The Green’s functions in the r.h.s. of eqs. (6.1b) and (6.1c) are calculated in the FDH 
scheme by setting ris = 4 as dictated^® and subtracting both poles and universal constants 
through the MS{x} operation. 

^®Here and in the following we explicitly write the values of Ua to highlight if contributions are regulated 
in FDH (ua = 4) or DReg (ua = n). In a fully general FDH scheme ria is an arbitrary parameter, but we 
restrict to the = 4 case. 
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If FDH amounts to a different renormalization scheme to MS, then the £-loop-CTs 
should satisfy the renormalizability condition (3.4), i.e. they should be proportional to 
the Born and local. We can proceed to calculate these at each order in a similar way 
to section 4, making use of information from the FDR vacuum part. At one loop, the 
calculation is simple and provides the exact same as FDR, indicating that FDH exhibits 
the same coupling shift from MS. At two loops the computation is more complicated but in 
general we can write the following expression for what should be the two-loop counterterms: 


(2-loop-CTs) = MS{GV+SV-EEI+CT}\n,=4 
- [GV+ SV-EEI+ GT] \n,=n- 


( 6 . 2 ) 


We now split SV and EEI as in eqs. (4.21) and (4.26). As Ug = n in the second line, we 
can make use of eq. (4.27). A direct calculation of the first line gives instead a striking 
different result for Ug = 4: the GT part cancels the b piece of the sub-vacuum, but not the 
EEI. In summary 


(2-loop-CTs) = MS{GV+SV^-EEIa-EEh}\n ,=4 

-[GV+SVl,-EEIa]\n^=n. 


(6.3) 


This is problematic for FDH as the EEI}, term does not in general satisfy the renormal¬ 
izability condition. That is, we see that for EEI}, / 0, FDH cannot be interpreted as a 
perturbative description of QFT. Nevertheless, as we have seen previously, the EEIj, con¬ 
tribution is zero for external gauge states and so the only expected problems come with 
external fermions, just as experienced in [7]. 

The form of equation (6.3) naturally suggests a new definition of FDH which sat¬ 
isfies the renormalizability condition. Let us consider changing the bare two-loop FDH 
correlation function in the following way^^ 


Gi 


bare I4 


G 


bare I 


=4 + EEI},\n^=4 


(6.4) 


This now allows us to write down a working FDH analogue of eq. (4.28) 

Z2,FDH'qO ^ m{GV+SVyEEIa}\ns=4 . , 

~[GV+SV;-EEIa]\ns=n. 

Here we call this modified definition FDH^ Its renormalization constants given 

in appendix A and correspond to the case (A, 5) = (0,0). From these we are able to 
calculate the analogous version of eq. (5.1) in this scheme, indicating that the QCD S-T 
identity is respected, even with external fermions^^. What’s more, we verify that in all cases 
the coupling constant shift agrees with the literature value [11], i.e. FDH' is equivalent 
to FDH when the latter scheme makes sense and provides consistent predictions in all the 
other cases. Together these results suggest that, at least off-shell, this definition does not 

is important to understand that the change we make in equation (6.4) can be read directly from the 
diagrams without the computation of counterterms. 

^®This is the important new case as the EEIs are zero for the correlation functions with external gauge 
states. 
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face the renormalization difficulties of the original FDH formulation. Furthermore, as we 
newly have control over the fermion sector, we are able to calculate a mass renormalization 
constant, and thereby a mass shift between MS and FDH' 


m; 


FDH' 


= I-Cf 


a 


MS 


m, 


MS 


dvr 


+ Cf 


a 


MS 


dvr 


12 2 d-' 


( 6 . 6 ) 


Our new definition offers a different perspective than the recently proposed approaches 
to the unitarity based difficulties of FDH. Roughly speaking, any solution requires a con¬ 
nection between internal and external states. In the approach of [12], evanescent operators 
are introduced, as in dimensional reduction, to make the external states behave like the 
internal ones, thereby introducing a series of evanescent couplings. Conversely, our ap¬ 
proach based upon sub-integration consistency makes internal states behave like external 
ones, and so does not require the introduction of any new fields whilst maintaining = d 
spin degrees of freedom. 


7. Summary and outlook 

As the FDR UV subtraction is consistently encoded in the definition of a four-dimensional 
and finite loop integration, the FDR approach to QFT does not require the introduc¬ 
tion of counterterms in the Lagrangian. In particular, an order-by-order renormalization 
is avoided: the perturbative order is computed by only looking at £-loop Feynman 
diagrams. 

We have proven, up to two loops in QCD, that FDR is equivalent to a particular 
renormalization scheme of Dimensional Regularization, dubbed Djfeg'', whose renormal¬ 
ization constants, extracted from two- and three-point vertices, obey the Slavnov-Taylor 
identities. Djfeg'' and MS are related by shifts in ag and niq that we have explicitly com¬ 
puted. The transition rules derived in this paper can be used to translate calculations of 
IR finite quantities from FDR to MS. 

During our analysis, we have identified a bottom-up interpretation of the failure of the 
naive FDH scheme beyond one loop. FDR provides a natural fix to this: it automatically 
generates the finite pieces needed to restore the renormalizability of FDH. These extra terms 
are computed without introducing e-scalars nor evanescent quantities in the Lagrangian 
and can be directly read off from two-loop diagrams. Including such FDR inspired terms 
in FDH defines a consistent renormalizable scheme that we have called FDH'. 

In this paper we have concentrated our focus on IR finite quantities. The possibility 
of consistently using FDR to regulate final state soft/collinear divergences at one-loop has 
been proven in ref. [15]. The study of FDR in the presence of NNLO IR singularities is left 
for future investigations. 
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A. The renormalization constants of the QCD vertices 


In the following we list the renormalization constants of the QCD correlation functions of 
figure 1 up to two loops in the Feynman-’t Hooft gauge 

Z,{\5) = 1 + (2) Z}{\,5) + (g)'z2(A,5). (A.l) 

The three renormalization schemes we are interested in are parametrized by the values of 
A and 6 given in table 1. Nc is the number of colours, Nf the number of active fermions, 

j\^2 _^ 

and as = as (A, <5) is the QCD coupling constant of each scheme. The constant 


/ = ^(Li2(e*f)-Li2(e-t)) 


-1.17195361... 


originates from the evaluation of the two-loop global vacuum in DReg, 


(A.2) 



JA FDR 
-L-'Reg 

MS 

FDH' 

A 

0 

1 

0 

(5 

1 

0 

0 


Table 1: A and 6 in Dneg'^, MS and FDHA 


Zgg = Zi{X,S) 
= 1 + 




dvr 

(S)^ 


ri 

[5 2 1 

17 

1 CO 

1 

1 CO 


+ - A) 


-—N^ + -NcNf 
12 6 ^ 


f^N^ + lxN^ + Nf^-C,--N, 


1 

-h- 
e 

+(l-A)^A2 + 5(liV2^IiV2/ 

+Nf (^-(1 - A)^^. + ,5 Qa, -^CF + f(^ 

Zee = Z2{X,S) 


--Cf - -Ne 
3 3 


= ' + ?- 

dvr e 




(^Y J _ 


-Ae 
2 

2 r 1 


V dvr 7 


-Ai + -NeNf 


1 

+- 

€ 


^a2 + - ^NeNf 

d8 d 2d ■' 


«(fivi+ InY +w/ - Yj 


(A.3) 


(A.4) 


- 16 ~ 



Zxpq, = Z3(A, 6) 


+ 


\ 471 


1 


NcCf + \cl 


+- 


1 


1 


17 


N,CF + \C^p + -NfCF 


+(1 — A) (^^NcCf — 2^f'^ + ^ ^^cCf + / 

+Nf (^-(1 - X)^Cf + 5(cf + ^CfJ 


-_CP-4N,Cf 


(A.5) 


Zggg = Z4{X,6) 

= 1 + ? 
4f 

4f ) 


+ 

1 

+- 

e 


ri 

r2 2 1 

- 

CO 1 

1 

CO 1 

U 


1 r 13 ^. 

a2 [ 8 " 


+ 3^c(l — A) 


r^^ciV/ 


25 


|jV|+lAJVf + JV,^-Cf-^iV, 


+ (1 - A)^ivf + « + jA'"/ 

+JV^ 1_(1 _ a)Cf + S - Icf + / 


_lCir--lVc 
3 6 



(A.6) 


Zgcc = ^5(A, 5) 
as 1 


= 1 + 


47r e 


--iVc 

2 


+(t^) W 


as 


4f 


-Ni 


1 

+ - 
e 




+ 5 


3 o 3 n 
--Nl - -Nlf 

8 c 4 CJ 


(A.7) 


Zg'^’^ — Zq{X, 5) 
as 


= 1 + 

+ 


1 


47r 1^ e 
/as\2 f 1 
V 4f 


[-Cf - iVc] + ^^.(1 - A) 


25 n 1 n 1 
-1 + + -Cl - -N,Nf 


1 

+- 

e 


f - - ^A^c + A f-i + 


2 4 

1 1 
2 + 4- 


+ A'/lA;Vo + icF 




+A^/ (-(1 - ^)Icf + 5(cf- ^iV, + / QiV, + ^Cf 



(A.8) 
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The above results can be used to derive the renormalization constant of the QCD coupling 
through eqs. (5.2) 




-—Nc + -Nf 

3 3 •' 


+ 




Vdvr/ 1 


^^0 + >1 - TJV.JV/ 


- - A) 

44 


9 


9 


9 


1 

+ - 

e 


-—N^ - —\Nl + Nf(—Nc + Cf + -\N^ 
9 c 9 c M 9 ^ -^9 


+ (-(1 - A)^iV2 + 5 ) 

+Nf ((1 - \)Cf + 5 Qcf -N, + f QiV, + Icf) ) 


(A.9) 


Finally, computing G\,^, with niq ^ 0 gives the renormalization constant associated with 
the quark mass 

Zrriq = 1 + /— [—SCi?] + Ci?(l — A) 


+ 

1 

+- 

e 




2 r 1 


dvr/ 1 
15 


yiVcCp + ^-Cl - NjCf 
85 


-^Cl - ^N,Cf + A [3Cl - N,Cf) + ^NfCF 
+(1 — A) -^^cCf'^ — S ^—Cp + —NcCf + / ^—Cp + 9NcCf 

+Nf (-(1 - X)\cf + 4^^/) }• (A-10) 

B. A two-loop diagram 

As an example of the algorithm we use to extract and compute the FDR vacuum we work 
out in detail one among the two-loop diagrams contributing to in figure 1, namely the 
ghost-loop correction depicted in figure 4. 

The integrand of the corresponding am¬ 
plitude is proportional to 




(^1 • qi2){qi + ki) ■ q2 


qtDiqlqj^ 


(B.l) 



with qi 2 — qi + q2 and Di — {qi + ki)'^. To Figure 4: Two-loop diagram contributing to 
read off J{qi, q 2 , from J{qi, ( 72 ) we apply the ghost self-energy, 
the shift of eq. (2.2) in both the numerator 
and the denominator and simplify reducible 
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numerators. This gives 

2^ _ (fei • q2f kj (fci • g2) 1 (fei • g2) 1 1 (fci • ^2) 1 1 

91 , 92 ,^ _ 4 ^^_ 2 - 2 ^ 2 qfDiqlql^ ‘^^ 1 ^ 2^12 ^ QiDiq^ qfDiql 

ki 1 _ 1 (fci • 92 ) _ 1 1 fcf 1 _ (fci • 92 ) 

4 g4/)ig22 2 q\Diql2 ^iDiqlql2 9f^i9i9?2 

11 11 11 11 

4 91^2912 4 gf Z)ig| 4 gf i)igf2 4 qlDiq\.^' 

Notice that the integrand manipulations we have performed so far are allowed both in FDR 
and DReg. In the following we concentrate on the first tensor 

=glis' 

and explicitly derive the splitting 

Jo{qi,q 2 , = [«4o,INF(9l,<?2,li^)] + «4 o,f(91)92,M^) (B.4) 


m2 

(B.2) 


needed to define the FDR integral 

J[d^qi][d‘^q2] Jo(gi, ^2,1^^) = lim J (i^gi(i^g2 ^o,f(^i, 92, (B.5) 

The other terms in eq. (B.2) are easier and can be treated analogously. 

To analyze the UV behaviour of a two-loop integrand J we introduce the four operators 
(io(«4), di{J), d 2 {J) and di 2 {J)-, which indicate how the integral over it behaves for large 
values of the integration momenta when 91 —)■ 00 and 92 —)• 00 independently {do{J)) or 
when 9i is fixed {di{J)). Thus, UV divergences occur when di{J) > 0 for some i. In our 
case, as 


do{Jo) — di{Jo) — 2, d2{Jo) < 0, di2{Jo) < 0, (B.6) 


Jo(91,92,1^^) has a logarithmic global UV divergence and is quadratically divergent in one 
of its sub-integrations. We now apply twice the identity 


and rewrite 


with 



1 

-2 

1 

9i + 2(91 • 92) 

-2 -2 


(B.7) 


9i2 

92 

92912 



^2,J^) = 

[•fci(9i. 

92, J^)] 

— [•fc2 

(91,92, J^)] 

+ 'fc3(9l, 

92, (B.8) 

[•fci (91,92, 

J^)] = 

1 

9f^i 

9? 92 
. ^2 _ 



(B.9) 

[•fc2(9l,92. 

/^^)] = 

9? 

9^^i 

9? 92 

to 

biT 

92 92 92 

, (B.IO) 

•fc3(9l,92 

= 

9^92(9? +2(91 -92))^ 
9f^i^9f2 


(B.ll) 
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We see that 92 ;/^^)] and [^ 2 (<?i) 92 ;/^^)] are factorizable integrands in which the 

^ 2—^00 behaviour is fully parametrized in terms of divergent integrands depending only 
on fj?. Therefore, they belong to [Jo,inf(9i, 92 , • In addition, since 


'io(-l3) = 0, (ii(J3)=0, d2{J‘i)=Q, di2{Ji,)<Q, 

further infinities need to be extracted from it, that is achieved by rewriting 

1 _ 1 + 2{p- qi) 

W Wi ’ 

which gives 

>^3(91,92,^^) = [«A( 9 i, 92 ,M^)] - <^5(91,92,^^) - «/ 6 ( 9 i, 92 ,^^), 


(B.12) 


(B.13) 


(B.14) 


where 


[«/4(91,92,^^)] = 
- 15 ( 91 ,92,= 4 
-16(91,92,/!^) = 


9^92(9?+ 2(91 -92))' 


me,2 


9292(91-92) (p +2(p-91)) 
^^1^912 

92 92 (9t + 49 i(9i_- 92))(p^ + 2(p • 91 )) 
^-^192912 


(B.15) 


Thus, the non-factorizable integrand [- 14 ( 91 , 92 , contributes to [-Io,inf( 9 i, 92 , and 
Je{qi,q 2 , is UV convergent. Furthermore, since 


do{J5)<0, di{J5) — 0, d2{J5)<0, di2{J5)<0, (B.16) 

a logarithmic sub-divergence is still present in - 15 ( 91 ,92, that gets separated when 
applying once more eq. (B.7) 

- 15 ( 91 ,92, = [- 17 ( 91 ,92, - J8{qi,q2,P^), (B.17) 


with 




a P 1 S 

92 92 92 92 
92 


Mqi,Q2,P ) = 4 


21 _ ,9f92(9i-92)^(r + 2(p-9i))(9f-F2(gi-g2)) 


^-Di 92912 

In summary, Jo(9i, 92,should be split as follows: 
-Io,f( 9 i, 92,11^) = -13(91, 92,11^) - 


(B.18) 

(B.19) 


[-Io,inf( 9 i, 92,I1 ^)] = [-14(91,92,11^)] + [-Ii( 9 i, 92,M^)] - [-^2(91,92, Ii^)] - [-17(91,92, Ii^)] • 

In the rest of the appendix we analyze the four terms contributing to [-Io,inf(9i, 92, Ii^)] to 
establish their connection with GV 2 {p) and SV 2 {p) in eqs. (4.12)-(4.17). 
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It gives a contribution to GV 2 {n). By using integration-by-parts one finds, for any 
value of n, 


lim 

l-L —^0 ^ 


d^qi d^q2 r , 2 n 1 

^ [J4{qi^q2,d‘ )J 


dn d'K 


= 5 "^ lim 


M—, 


d^qi d»g2 


4 

3n 


-4 -4 


1 4 \ 

’ 6 ~ 


19192912 


(B.20) 


• [«/i(9i,92,M^)]; 

Since it is proportional to fj?, it vanishes in the limit /U^ —)• 0 we are interested in. 


• [Mqi,q2,d^)]- 

The the odd-rank tensor in the r.h.s. of eq. (B.IO) gives zero upon integration. On 
the other hand —)■ 0 is allowed in the qi integrand multiplying the rank-two tensor. 
Thus 


lim 

f-L —^0 


d”gi d'^q2 p ^ , 2 m 

-3^-327 [^ 2 ( 91 ,92,/i )J 

H'R Mr 


^0/3 

-lim 

4 At-i-O, 


d^qi d^q2 


1 

L92J 




(B.21) 


The first term in the curly bracket contributes to SV 2 {iJ.) with A '2 = — 2 I 2 = g°‘^ 1'^, 
while the second to GV 2 {g) with F 22 = 0 and F 21 = 


• [«/7(9i,92,M^)]: 

The integrand depending on qi in eq. (B.18) behaves as l/q\ when qi —)■ 0. To 
disentangle the term which develops a In g? we rewrite^® 


917915(4’^ + 2(p • gi)) _ 917915 917915 _ 917915 917915 


^T>i 


9i 


qlDi 




9t^i 


which results in 


[«/7(9i,92,M^)] = 4 

Since the integral over 91 is UV finite 

d"9i d'^q2 


( 

Ql'yQlS 

917915 \ 

1 - 

1_ 

V 

[ 91 J 

qfDiJ 

92 J 


I U mu H2 r , 2\1 

lim / ^ 27 ^^ [<47(91,92,/! )J 


/Xr^^r 


(B.22) 


(B.23) 


(B.24) 


only contributes to 5 V 2 (/x) • By using tensor reduction and integration-by-parts iden¬ 
tities one finds yl 2 = and 7 I 2 = ^9^'^ for this term. 

^®As mentioned in section 4, this separation is in general achieved by using the propagator identity in 
eq. (4.15). 
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C. Sub-prescription example 


Here we discuss a further example of the sub-prescription in order to aid the understanding 
of the reader for future FDR calculations. Consider the FDR integral 

7/3(^2 + Vl)7a(^i2 + + Vl)7“ 


/ 


[d\i][d%]- 



92 


Figure 5: Diagram contributing to Gg. 
fermion momenta follows the fermion line. 


qlqlDiD2Di2 

corresponding to the contribution to Gg in figure 5, where Di = qf + k'f + 2{qi ■ ki) 
We wish to discuss how to extract the 
EEI resulting from the sub-prescription, 
so as in figure 2 we consider the un-promoted 
numerator in order to find the relevant terms. 

In this diagram we have two sub-divergences, 
one for fixed qi and another fixed q 2 - The 
terms from the sub-prescription of each can 
be extracted considering the sub-divergences 
independently. First we shall consider qi 
hxed. Let us disconnect the divergent sub¬ 
diagram, as in figure 6. 

The numerator with its appropriate 
hatting reads 

N = 7/3(^2 + ^i)7«(^i 2 + ^i)7^(^i + ¥l)7" 

(C.2) 

and the numerator terms which give loga¬ 
rithmic divergences in q 2 are 
Ar(2) = 

= + Vi) - 4^2^^! + 

(C.3) 

The sub-prescription gives 

^(2) ^^(2) ^4^21 

while the global prescription requires no change in 

Ar(2) ^ jv{2)^ 


(C.l) 


The 



Figure 6: The same diagram of figure 5 from the 
point of view of the q 2 sub-integration. 


(C.4) 


(C.5) 


Thus, we find the resulting contribution by subtracting the (zero) global promotion and 
adding in the sub-promotion, leading to an EEI of the form^*^ 

AK^l + Vl) _^._2 ^^4.1 Vl 


EEI = A [d^qi][d%]: 


= ITT' 


[d^qi]-_ 


(C.6) 


' q\qlDiD2Di2 J q\Di 

When we move to perform the sub-prescription in the second sub-divergence, i.e. at fixed 
52 , we make a similar treatment and find that the sub-prescription gives an identical con¬ 
tribution due to the symmetry of the diagram. 

^°See appendix D. 
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D. Computing EEIs 


Given the important role played by the EEIs in the consistency of FDR we explicitly 
compute the extra-extra integrals in eqs. (4.36) and (C.6): 

4„ 1 + 


EEh= /[d"(?i][d"g2]=2=^ 


+kl + 2qi ■ ki){ql + kl + 2q2 ■ ki){qf2 + kj + 2qi2 ■ ki 


As a first step, we need the related one-loop extra (sub-)integrals 


hipl) = 




h= j [d'^q] 




(PD 1 D 2 ' 


(D.l) 


(D.2) 


where Di = {q + Pi)"^ and Di = Di — fj?. To calculate them we start from the FDR defining 
expansions of their integrands with —)■ q^ 
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q^Dl 
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'q^iq-piY 
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'q^{q-pif' 



. <f . 


q^ 


+ Ci(g), 


+ F2{q), 


(D.3) 


^DiD2 

where Fi{q) and ^ 2 ( 9 ) are UV convergent. Ii{pi) and I 2 are defined [6] as the difference 
between the l.h.s. of eq. (D.3) and the UV divergent part computed by changing back 


in the numerator 

IM) = limy [ rq 


I 2 = lim fj, / (Fq 

fi^O J 
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Pi , n (g-Pl) Ad-Plf 


q'^Df q'^ (f 

1 1 \ ivr^ 

q‘^DiD2 “ Fy ^ ^ 




■ 2 P 1 

= -T, 
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(D.4) 


Therefore we obtain 
EEh = 


EEl2 = 




= - 


+ 2gi • ki) 

+ Vi 


ITT 


12 




+ 2gi • /ci)’ 


l2 = 


ITT 




1 


qfiqj + kf + 2qi ■ ki)' 


(D.5) 


kl + 2qi ■ ki) 

Notice the replacement hiqf) —)■ hiof) in accordance with the global prescription. 

Finally, we point out the difference between hatting and not hatting p^\ 2 , i-e. the 
inequivalence between EEIs and two-loop extra integrals. For instance 

/. ki 


EEh = -- 


12 


ln^-2 

Pi 


while 


f.A ir,4 1 P\{iil 


p\{il±_tl)__ ^ _T_ 

k\ + 2qi ■ ki) 12 


/ kf 5 16 


In 


Pi 


+ o + -rrf • 


(D.6) 


(D.7) 


- 23 - 



References 


[1] N. Bogoliubov and O. Parasiuk, On the Multiplication of the causal function in the quantum 
theory of fields, Acta Math. 97 (1957) 227-266. 

[2] K. Hepp, Proof of the Bogolyubov-Parasiuk theorem on renormalization, Commun.Math.Phys. 
2 (1966) 301-326. 

[3] W. Zimmermann, Convergence of Bogolyubov’s method of renormalization in momentum 
space, Commun.Math.Phys. 15 (1969) 208-234. 

[4] G. ’t Hooft and M. Veltman, Regularization and Renormalization of Cauge Fields, 

Nucl.Phys. B44 (1972) 189-213. 

[5] J. C. Collins, Renormalization. Cambridge University Press, 1984. 

[6] R. Pittau, A four-dimensional approach to quantum field theories, JHEP 1211 (2012) 151, 
[arXiv; 1208.5457]. 

[7] W. B. Kilgore, Regularization Schemes and Higher Order Corrections, Phys.Rev. D83 (2011) 
114005, [arXiv: 1102.5353]. 

[8] R. Boughezal, K. Melnikov, and F. Petriello, The four-dimensional helicity scheme and 
dimensional reconstruction, Phys.Rev. D84 (2011) 034044, [arXiv: 1106.5520]. 

[9] W. B. Kilgore, The Four Dimensional Helicity Scheme Beyond One Loop, Phys.Rev. D86 
(2012) 014019, [arXiv: 1205.4015]. 

[10] Z. Bern and D. A. Kosower, The Computation of loop amplitudes in gauge theories, 
Nucl.Phys. B379 (1992) 451-561. 

[11] Z. Bern, A. De Freitas, L. J. Dixon, and H. Wong, Supersymmetric regularization, two loop 
QCD amplitudes and coupling shifts, Phys.Rev. D66 (2002) 085002, [hep-ph/0202271]. 

[12] A. Broggio, C. Gnendiger, A. Signer, D. Stckinger, and A. Visconti, Computation of H ^ gg 
in FDH and DRED: renormalization, operator mixing, and explicit two-loop results, 

arXiv:1503.0910. 

[13] A. Donati and R. Pittau, Gauge invariance at work in EDR: H —>■ gamma gamma, 

arXiv:1302.5668. 

[14] A. M. Donati and R. Pittau, FDR, an easier way to NNLO calculations: a two-loop case 
study, Eur.Phys.J. C74 (2014) 2864, [arXiv: 1311.3551]. 

[15] R. Pittau, QCD corrections to H ^ gg in FDR, Eur.Phys.J. C74 (2014) 2686, 

[arXiv: 1307.0705]. 

[16] R. Pittau, On the predictivity of the non-renormalizable quantum field theories, Portsch.Phys. 
63 (2015) 132-141, [arXiv: 1305.0419]. 

[17] R. Pittau, Integration-by-parts identities in FDR, arXiv : 1408.5345. 

[18] P. Nogueira, Automatic Feynman graph generation, J.Comput.Phys. 105 (1993) 279-289. 

[19] J. Kuipers, T. Ueda, J. Vermaseren, and J. Vollinga, FORM version 4-0, 
Comput.Phys.Commun. 184 (2013) 1453-1467, [arXiv: 1203.6543]. 

[20] E. Egorian and O. Tarasov, Two Loop Renormalization of the QCD in an Arbitrary Gauge, 
Teor.Mat.Fiz. 41 (1979) 26-32. 


- 24 - 



[21] T. Muta, Foundations of Quantum Chromodynamics: An Introduction to Perturbative 
Methods in Gauge Theories, (3rd ed.), . 

[22] L. N. Mihaila, J. Salomon, and M. Steinhauser, Renormalization constants and beta functions 
for the gauge couplings of the Standard Model to three-loop order, Phys.Rev. D86 (2012) 
096008, [arXiv; 1208.3357], 

[23] S. Larin and J. Vermaseren, The Three loop QCD Beta function and anomalous dimensions, 
Rhys.Lett. B303 (1993) 334-336, [hep-pli/9302208]. 

[24] J. Vermaseren, Axodraw, Comput.Phys.Commun. 83 (1994) 45-58. 


- 25 - 



